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Abstract — Sparse intersymbol-interference (ISI) channels are 
encountered in a variety of high-data-rate communication sys- 
tems. Such channels have a large channel memory length, but only 
a small number of significant channel coefficients. In this paper, 
trellis-based equalization of sparse ISI channels is revisited. Due 
to the large channel memory length, the complexity of maximum- 
likelihood detection, e.g., by means of the Viterbi algorithm (VA), 
is normally prohibitive. In the first part of the paper, a unified 
framework based on factor graphs is presented for complexity re- 
duction without loss of optimality. In this new context, two known 
reduced-complexity algorithms for sparse ISI channels are reca- 
pitulated: The multi-treUis VA (M-VA) and the paraUel-trellis VA 
(P-VA). It is shown that the M-VA, although claimed, does not lead 
to a reduced computational complexity. The P-VA, on the other 
hand, leads to a significant complexity reduction, but can only be 
applied for a certain class of sparse channels. In the second part 
of the paper, a unified approach is investigated to tackle general 
sparse channels: It is shown that the use of a linear filter at the 
' receiver renders the application of standard reduced-state trellis- 
' based equalizer algorithms feasible, without significant loss of op- 
timality. Numerical results verify the efficiency of the proposed 
^ receiver structure. 

I Index Terms — Trellis-based equalization, sparse ISI channels, 
I complexity reduction, prefiltering. 



I. Introduction 

SPARSE intersymbol-interference (ISI) channels are en- 
countered in a wide range of communication systems, such 
as high-data-rate mobile radio systems (especially in hilly ter- 
rain), wireline systems, or aeronautical/ satellite communica- 
tion systems. For mobile radio applications, fading channels 
are of particular interest. The equivalent discrete-time chan- 
nel impulse response (CIR) of a sparse ISI channel has a large 
channel memory length L, but only a small number of signifi- 
cant channel coefficients. 

Due to the large memory length, equalization of sparse ISI 
channels is a demanding task. The topics of linear equalization 
and decision-feedback equalization for sparse ISI channels are, 
e.g., addressed in |1 1, where the sparse structure of the chan- 
nel is explicitly utilized for the design of the corresponding 
finite-impulse-response (FIR) filter(s). Trellis-based equaliza- 
tion for sparse channels is addressed in |2 1-|4|. The complexity 
in terms of trellis states of an optimal trellis-based equalizer, 
based on the Viterbi algorithm (VA) ||5] or the Bahl-Cocke- 
Jelinek-Raviv algorithm (BCJRA) |6^, is normally prohibitive 
for sparse ISI channels, because it grows exponentially with 
the channel memory length L^. However, reduced-complexity 

'^The VA is optimal in tlie sense of maximum-likeliliood sequence estimation 
(MLSE) and tlie BCJRA in tlie sense of maximum a-posteriori (MAP) symbol- 
by-symbol estimation. Both algorithms operate on the same trellis diagram. All 
statements concerning complexity hold for both the VA and the BCJRA. 



algorithms can be derived by exploiting the sparseness of the 
channel. 

In ||2|, it is observed that given a sparse channel, there is 
only a comparably small number of possible branch metrics 
within each trellis segment. By avoiding to compute the same 
branch metric several times, the computational complexity is 
reduced significantly without loss of optimality. However, the 
complexity in terms of trellis states remains the same and thus 
the storage expense. As an alternative, another equalizer con- 
cept coined multi-trellis Viterbi algorithm (M-VA) is proposed 
in |2| that is based on multiple parallel irregular trellises (i.e., 
time-variant trellises). The M-VA is claimed to have a signifi- 
cantly reduced computational complexity and storage expense 
without (much) loss of optimality. 

A particularly simple solution to reduce the complexity of 
the conventional VA without loss of optimality can be found 
in 1 3 1: The parallel-trellis Viterbi algorithm (P-VA) is based on 
multiple parallel regular trellises. However, an application of 
the P-VA is only possible for a certain class of sparse channels 
having a so-called zero-pad structure. In order to tackle more 
general sparse channels with a CIR close to a zero-pad channel, 
it is proposed in |3| to exchange tentative decisions between 
the parallel trellises and thus cancel residual ISI. This modified 
version of the P-VA is, however, suboptimal and is denoted as 
sub -P-VA in the sequel. 

A generalization of the P-VA and the sub-P-VA can be found 
in 1 4 1, where corresponding algorithms based on the BCJRA 
are presented. These are in the sequel denoted as parallel-trellis 
BCJR algorithms (P-BCJRA and sub-P-BCJRA, respectively). 
Some interesting enhancements of the (sub-)P-BCJRA are also 
discussed in |4|. Specifically, it is shown that the performance 
of the sub-P-BCJRA can be improved by means of minimum- 
phase prefiltering |7|, |8| at the receiver. A specific FIR ap- 
proximation of the infinite-length linear minimum-phase filter 
is used, which preserves the sparse structure of the channel. 
This guarantees that the sub-P-BCJRA can still be applied after 
the prefiltering. 

Alternatives to trellis-based equalization are the tree-based 
LISS algorithm |9| and the Joint Gaussian (JG) approach 
in 1 10 1 . In this paper, trellis-based equalization for sparse ISI 
channels is revisited. A unified framework based on factor 
graphs fln is presented in SectionlUlfor complexity reduction 
without loss of optimality, and the M-VA as well as the P-VA 
are recapitulated in this new context. Specifically, it is shown 
that the M-VA does, in fact, not lead to a reduction of compu- 
tational complexity, compared to the conventional VA. More- 
over it is illustrated, why the optimal P-VA can only be applied 
for zero-pad channels. In order to equalize general sparse ISI 
channels, a simple alternative to the sub-P-VA/ sub-P-BCJRA 



is investigated in Section Hill For this purpose, the idea in |4| 
to employ prefiltering at the receiver is picked up. It is demon- 
strated that the use of a Unear minimum-phase filter renders the 
application of reduced-state equalizers such as |12|, |13| fea- 
sible, without significant loss of optimality. The proposed re- 
ceiver structure is notably simple: The employed equalizer al- 
gorithms are standard, i.e., not specifically designed for sparse 
channel, because the sparse channel structure is normally lost 
after prefiltering. Solely the linear filter is adjusted to the cur- 
rent CIR (which is particularly favorable with regard to fading 
channels), where the filter coefficients can be computed accord- 
ing to standard techniques available in the literature |7|, |8|. In 
order to illustrate the efficiency of the proposed receiver struc- 
ture, numerical results are presented for various types of sparse 
ISI channels. Bit error rates (BERs) are achieved that deviate 
only 1-2 dB from the matched filter bound (at a BER of 10^"^). 
To the authors' best knowledge, similar performance studies for 
prefiltering in the case of sparse ISI channels have not yet been 
presented in the literature. 

II. Complexity Reduction without Loss of 
Optimality 

A general sparse ISI channel has a comparably large chan- 
nel memory length L, but only a small number of significant 
channel coefficients hg, g = 0, G L, according to 

h := [ho hi ... hcf, (1) 

/o zeros fi zeros /g - 1 zeros 

with fi>0 integer for all i and L = Yl'i=o{fi + !)■ ^ sparse 
ISI channel, for which /o = /i = ■■■ = fo-i —'■ f holds, is 
referred to as zero-pad channel OJ. 

Throughout this paper, complex baseband notation is used. 
The fc-th transmitted A/-ary data symbol is denoted as x[k], 
where k is the time index. A hypothesis for x[k] is denoted 
by x[k] and a hard decision by x[k]. For simplicity, the channel 
coefficients are assumed to be constant over an entire block of 
data symbols (block length N > L). The equivalent discrete- 
time channel model is given by 

G 

y[k] = hQx[k] + ^ hgx[k — dg] + n[k] , (2) 

9=1 

where y[k] denotes the fc-th received sample and n[k] the fc-th 
sample of a complex additive white Gaussian noise (AWGN) 
process with zero mean and variance cr^. Moreover, 

a 

dg ■■= + l<g<G, (3) 

1=1 

denotes the position of hg in h. 

In the sequel, the channel vector h is assumed to be known 
at the receiver. Moreover, an M -ary alphabet for the data sym- 
bols is assumed. The complexity in terms of trellis states of 
the conventional Viterbi/ BCJR algorithm is given by 0{M^} 
and is thus normally prohibitive. Given a zero-pad channel, it is 
proposed in 1 3 1 to decompose the conventional trellis diagram 
with = A/(-''+i)'^ states into (/ + !) parallel regular trel- 
lises, each having only Af*^ states. As will be shown in the 



sequel, a decomposition into multiple parallel regular trellises 
is not possible in the case of a more general sparse ISI channel. 

In order to decompose a given trellis diagram into multi- 
ple parallel trellises, the following question is of central inter- 
est: Which symbol decisions x[k] are influenced by a certain 
symbol hypothesis i[fco]? Suppose, a certain decision x[ki] is 
not influenced by the hypothesis .T[fco], where fco and fci are 
two arbitrary (but fixed) time indices. Furthermore, let the set 
Xkf, :^{x[k] \x[k] depends on x[fco]} contain all decisions in- 
fluenced by x[ko] and the set Xki all decisions influenced by 
x[ki]. If these two sets are disjoint, i.e. Xkg n Xk^ =0, the hy- 
potheses x[ko] and x [ki] can be accommodated in separate trel- 
lis diagrams without loss of optimality. In this case, a decom- 
position of the overall trellis diagram into (at least two) parallel 
regular trellises is possible. 

This fact is illustrated in the sequel for two example channels 
(L = 8 and G = 2 in both cases). The first channel is charac- 
terized by a CIR h(i) := [/lo /ii /la]"^ and 
the second by a CIR h^^) [ /iq /ii /i2 
For the time being, a particular symbol hypothesis x[ko] is con- 
sidered, and it is assumed that hard decisions x[k] are already 
available for all time indices k < ko. Moreover, it is assumed 
that the hypothesis x[ko] does not have any impact on decisions 
X [k] with fc > fco + DL {D > integer; in the example, we 
consider the case D ~ 2)^. Fig. 1 shows a diagram for the 
first example CIR. The diagram may be interpreted as a factor 
graph 1 1 1 1 and illustrates the dependencies between hypothe- 
sis i[fco] and all decisions x[k], ko < k < ko + DL. It can be 
seen from (|3 that only the received samples j/[fco], j/[fco + 6], 
and j/[fco + 8] are directly influenced by the data symbol x[ko]. 
Therefore, there is a direct dependency between the hypothesis 
x[ko] and the decisions a;[fco], i[fco + 6], and a;[fco + 8]. The re- 
ceived sample y[fco+8], for example, is also directly influenced 
by the data symbol x[fco + 2]. Correspondingly, there is as well 
a dependency between x[ko] and the decision 2;[fco + 2]. The 
data symbols a;[fco+6] and a;[fco+8] again directly influence the 
received samples y[ko+l2], j/[fco+14], and y[fco+16], and so on. 

Including all dependencies, one obtains the second graph of 
Fig. 1. As can be seen, there is a dependency between a; [fco] 
and all decisions x[ko + 2i'], where 1^ = 0,1, [DL/2\. Con- 
sequently, in this example it is possible to decompose the con- 
ventional trellis diagram into two parallel regular trellises, one 
comprising only the time indices ko+2v and the other one com- 
prising the time indices fcp + 1 + 2v. While the conventional 
trellis diagram has trellis states, there are only A/"* states 
in each of the two parallel trellises. (Moreover, a single trellis 
segment in the parallel trellises spans two consecutive time in- 
dices.) This result is in accordance with |3 1, since the CIR h^^^ 
constitutes a zero-pad channel [h'o h'l h'2 /13 h'/^\^ , 
where G = 4, / = 1, and h'l = h^ ~ 0. Generally spoken, a 
decomposition of a conventional trellis diagram into multiple 
parallel regular trellises is possible, if all non-zero channel co- 
efficients of the sparse ISI channel are on a zero-pad grid with 
/ > 1. In this case, the optimal P-VA can be applied^; otherwise 
one has to resort to the sub-P-VA or to alternative solutions such 
as the M-VA. 

^This con'esponds to the assumption that a VA with a decision delay of DL 
symbol durations is optimal in the sense of MLSE. 

■^The P-VA is still optimal in the sense of MLSE. The finite decision delay 
DL is not required and has only been introduced for illustrative purposes. 
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Fig. 1 

Dependencies BETWEEN SYMBOL HYPOTHESIS x[feo] AND SUBSEQUENT DECISIONS Example CIR h(^) ;= [/lo hi ?t2 J""" 

(L = 8, G = 2). 



The second CIR h*^^^ constitutes a counter example. Here, 
the symbol hypothesis x[ko] influences all decisions x[k'], 
ko<k' < kQ + DL (not depicted due to space limitation). Con- 
sequently, a decomposition of the conventional trellis diagram 
into multiple parallel regular trellises is not possible here. Still, 
a decomposition into multiple parallel irregular trellises is pos- 
sible, as proposed in 1 2 1 in the context of the M-VA. However, 
a significant subset of the dependencies resulting from the cor- 
responding factor graph is neglected in |2|. If all dependencies 
are taken into account, the M-VA does not yield any complexity 
advantage over the conventional VA. The computational com- 
plexities in terms of the overall number of branch metrics that 
have to be computed for a single decision x[ko] are stated in Ta- 
ble 1, for the conventional VA, the P-VA (example CIR h^^)), 
and the M-VA (example CIR h'^)). 

III. Prefiltering for Sparse Channels 

The preceding section has shown that trellis-based equaliza- 
tion of a general sparse ISI channel is quite a demanding task: 
An application of the optimal P-VA (or the P-BCJRA) is only 
possible for zero-pad channels. In the case of a more general 
sparse channel, the suboptimal sub-P-VA with residual ISI can- 
cellation can be used. However, for a good performance the 
CIR should at least be close to a zero-pad structure (Sj. The 
M-VA, on the other hand, was designed for sparse channels 
with a general structure, but does not offer any complexity ad- 
vantage over the conventional VA if all dependencies between 
the individual symbol hypotheses are taken into account. 

In order to tackle general sparse ISI channels, a simple alter- 
native to the sub-P-VA/ sub-P-BCJRA is proposed in the sequel: 
It is demonstrated that the use of a linear minimum-phase filter 
at the receiver renders the application of standard reduced-state 
equalizer algorithms feasible, without significant loss of opti- 
mality. The receiver structure under consideration is illustrated 



TABLE I 

Computational complexities: Viterbi algorithm (VA), 

PARALLEL-TRELLIS VA (P-VA), AND MULTI-TRELLIS VA (M-VA). 



Conventional VA, 

any CIR withL = 8 


P-VA, 

example CIR h'^' 


M-VA, 

example CIR h'^' 




0{2M'^} 


C»{3 M"* + 2 M« 
-h2M'^ + 2M® 
+ 2M^ + 2M^ 
+ 2M^ + M'^} 



in Fig. 2, where z [k] denotes the fc-th received sample after pre- 
filtering and hmin the filtered CIR. Within the scope of this pa- 
per, the ideal linear minimum-phase filter is approximated by 
an FIR filter of length L-p, where the approach in |8| is used 
to calculate the filter coefficients. The resulting FIR filter ap- 
proximates a discrete-time whitened matched filter (WMF), i.e., 
the effect of noise coloring is negligible. The computational 
complexity of calculating the filter coefficients is 0{Lf L^), 
i.e. only linear with respect to the filter length. Therefore, 
comparably large filter lengths are feasible. With regard to the 
trellis-based equalizer, we focus on delayed decision-feedback 
sequence estimation (DDFSE) 1 13| in the sequel. The number 
of trellis states in the DDFSE equalizer is , where K<^L is 
a design parameter In order to obtain a complexity that is sim- 
ilar to that of the (sub-)P-VA/ P-BCJRA equalizer, one should 
choose K such that^ 

K < logM(/ + l) + G. (4) 

In the following section, it is shown that the sparse channel 
structure is normally lost after prefiltering. Afterwards, numer- 

"^In order to find an appropriate value for K in the case of a general sparse 
ISI channel, one has to find an underlying zero-pad channel with a structure as 
close as possible to the CIR under consideration. 
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Fig. 2 

Receiver structure under consideration. 



ical results are presented for various examples to demonstrate 
the efficiency of the proposed receiver structure. 

A. Channel Structure After Prefiltering 

The sparse structure of a given CIR h is normally lost after 
minimum-phase prefiltering. For example, consider the CIR 
h= [Vol Vol VOA]^ (i.e. ||h||2 = 1), and let i?(z) 
denote the z-transform of h. The zeros of H{z) are given by 

• zo 1 2 = 0.69 ± j0.80 (\zo 1 2I = 1.06) and 

. ^0 3,4 = -0.69 ± jO.56 (|zo3,4| = 0.89). 
Correspondingly, h is neither minimum-phase nor maximum- 
phase^. The z-transform H^in{z) of the filtered CIR hnun is 
obtained as follows: Those zeros of H{z) that are inside or on 
the unit circle are retained for i7„iin(z), whereas the outside 
zeros are reflected into the unit circle. The resulting minimum- 
phase CIR is h^in = [0.79 0.12 - 0.02 0.20 0.56]'^, which 
is not sparse anymore. 

As opposed to this, in the special case of a zero-pad channel 
the sparse structure is always preserved after minimum-phase 
prefiltering: Let h := [ho hi ... hc]'^ denote a (non-sparse) 
CIR with z-transform H{z), and let hzp denote the correspond- 
ing CIR with memory length (/+ 1)G and z-transform Hzp (z), 
which results from inserting / zeros in between the individual 
coefficients of h. Furthermore, let zq, 1, zq, g denote the 
zeros of H{z). An insertion of / zeros in the time domain cor- 
responds to a transform z ^ ^^/(^"'"^^ in the z-domain, i.e., 
Hzp{z)^ H{zf+'^). This means, the (/+l)Gzeros of iJzp(2) 
are given by the (/+1) complex roots of zq, 1, zq, respec- 
tively. Consider a certain zero zo.g := tq, g exp(j c/?o, g) of 
H{z) that is outside the unit circle (r^^g > 1). This zero will 
lead to (/ + 1) zeros 



l/(/+l) 



'0,g 



exp J 



2ttX 



(5) 



of ffzp(^) (A = 0, 



./ + 1 

, /) that are located on a circle of ra 



dius rQ^g'^^^^ > 1, i.e., also outside the unit circle. By means 
of minimum-phase prefiltering, these zeros are reflected into 
the unit circle, i.e., the corresponding zeros of -ffzp,min(2) are 

given by 1/z^^* . Therefore, the sparse CIR structure is re- 
tained after minimum-phase prefiltering (with the same zero- 
pad grid), since the zeros of -ffzp,min(2:) are the (/ + 1) roots 
of the zeros of i/minl^). Specifically, the non-zero coefficients 
of hzp,min are given by the CIR h,nin. If the zeros of H{z) (or 
equivalently of i?zp,min(2)) are not too close to the unit circle, 
hmin is characterized by a significant energy concentration in 
the first channel coefficients. In this case the effective channel 
memory length of hzp is significantly reduced by minimum- 
phase prefiltering, namely by some multiples of (/ + 1), cf. Q. 

^Given a minimum-phase (maximum-phase) CIR h, all zeros of H{z) are 
inside (outside) the unit circle, i.e., |2o,i| !S 1 (1^0, i I ^ 1) for all *. 
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Fig. 3 

BER PERFORMANCE OF THE PROPOSED RECEIVER STRUCTURE IN THE 
CASE OF A STATIC SPARSE ISI CHANNEL. 



B. Numerical Results 



In the sequel, numerical results obtained by Monte-Carlo 
simulations are presented. To start with, a static sparse ISI 
channel is considered, and the BER performance of the pro- 
posed receiver structure is compared with the sub-P-BCJRA 
equalizer 12 • As an example, we consider here the CIR 
h = [ /lo . . . /i4 . . . /i7 . . . /ii5 with /iq = 0.87 and 
h4 — hj ^ /ii5 = 0.29 from 01, which has a general sparse 
structure (i.e., no zero-pad structure). The BER performance 
(binary antipodal transmission, M = 2) of the sub-P-BCJRA 
equalizer and the DDFSE equalizer with WMF is displayed in 
Fig. 3, as a function of E-^/Nq in dB, where £'b denotes the 
average energy per bit and Nq the single-sided noise power 
density {E\^/No l/Cn)- Due to the given channel memory 
length, the complexity of MLSE detection is prohibitive. As a 
reference curve, however, the matched filter bound (MFB) is in- 
cluded, which constitutes a lower bound on the BER of MLSE 
detection. The filter length of the WMF has been chosen as 
ip — 40. Since the channel is static, the filter has to be com- 
puted only once. When the parameter K is chosen as K — 4, 
the overall receiver complexity is approximately the same as for 
the sub-P-BCJRA equalizer. As can be seen in Fig. 3, the BER 
performance achieved by the proposed receiver structure is also 
comparable to the sub-P-BCJRA equaHzer At a BER of 10~^, 
the loss with respect to the MFB is only about 1 dB. At the ex- 
pense of a small loss due to residual ISI (0.5 dB at the same 
BER), the complexity of the DDFSE equalizer can be further 
reduced to K — 3. 

Next, we consider the case of a sparse Rayleigh fading chan- 
nel model, i.e., the channel coefficients hg (g ~0, ...,G) in ^ 
are now zero-mean complex Gaussian random variables with 
variance E{\hg\^} —: . It is assumed in the sequel that 
the individual channel coefficients are statistically independent. 
Moreover, block fading is considered for simplicity. As an ex- 
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Fig. 4 

BER PERFORMANCE OF THE PROPOSED RECEIVER STRUCTURE IN THE 
CASE OF A SPARSE RAYLEIGH FADING CHANNEL. 



ample, we assume a CIR with G = 3 and a power profile 

p := [<o 0^ <i (7^ 2 <3]T. (6) 

/ zeros 

Note that this CIR again does not have a zero-pad structure. By 
choosing different values for the parameter /, different chan- 
nel memory lengths L — f + 6 can be studied. To start with, 
consider a power profile with equal variances crj^ q = ... = 

cr^ 3 = 0.25. Fig. 4 shows the BER performance of the pro- 
posed receiver structure for binary transmission and three dif- 
ferent channel memory lengths L (solid lines: L = 6, dashed 
lines: L = 12, dotted Unes: L = 20). The filter length has been 
chosen as Lp = 20 (L = 6), Lp = 36 (L = 12), and Lp = 60 
{L = 20). As reference curves, the BER for flat Rayieigh fading 
(L = 0) is included as well as the MFB 1 14 Ch. 14.5]. In the 
case L = 6, MLSE detection is still feasible. As can be seen in 
Fig. 4, its performance is very close to the MFB. The DDFSE 
equalizer with if = 5 in conjunction with the WMF achieves 
a BER performance very close to MLSE detection (the loss at 
a BER of 10"'^ is only about 0.6 dB). Even when the chan- 
nel memory length is increased to L = 20, the BER curve of the 
DDFSE equalizer with WMF deviates only 2 dB from the MFB. 
Specifically, a significant gain compared to flat Rayieigh fading 
is achieved, i.e., a good portion of the inherent diversity (due 
to the independently fading channel coefficients) is captured. 
When the DDFSE equalizer is used without WMF, a significant 
performance loss occurs already for L — 6. For the larger chan- 
nel memory lengths, the influence of the WMF makes a dra- 
matic difference: The BER increases by several decades when 
the WMF is not used. 

Similar performance results have also been obtained for un- 
equal variances cr^ g. When the power profile p of the original 
CIR already exhibits an energy concentration in the first channel 
coefficients, the benefit of the WMF is smaller, but still signif- 
icant. As an alternative to the proposed receiver structure, we 
have also studied the use of a linear channel shortening filter 



(CSF) 1 15 1, which transforms a given CIR with memory length 
L into a CIR with a reduced memory length Ls< L. When Lg 
is chosen small enough, the application of the conventional VA 
is again feasible, operating on the shortened CIR. However, the 
performance of this receiver structure turned out to be inferior 
to the DDFSE equalizer with WMF 1 16 1. 

The concept of minimum-phase prefiltering for sparse ISI 
channels is also beneficial when using a tree-based equalization 
algorithm, such as the LISS algorithm |9|. In order to obtain a 
small complexity, the metrics of two competing paths that devi- 
ate closely to the root of the tree should differ as much as possi- 
ble. This is achieved by means of minimum-phase prefiltering, 
due to the energy concentration in the first channel coefficients 
of the filtered CIR. 

IV. Conclusions 

In this paper, trellis-based equalization of sparse inter- 
symbol-interference channels has been revisited. Due to the 
large memory length of sparse channels, efficient equalization 
with an acceptable complexity-performance trade-off is a de- 
manding task. With regard to complexity reduction, it has been 
demonstrated in which cases a decomposition of the conven- 
tional trellis diagram into multiple parallel trellises is possible 
without loss of optimality. In order to tackle general sparse 
channels, a receiver structure with a linear filter and a reduced- 
complexity trellis-based equalizer has been studied. The em- 
ployed equalizer algorithm is standard, i.e., not specifically de- 
signed for sparse channels, because the sparse channel structure 
is normally lost after prefiltering. The coefficients of the linear 
filter can be computed using standard techniques from the liter- 
ature. By means of numerical results, the efficiency of the pro- 
posed receiver structure has been demonstrated, both for static 
and fading channels. 
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Abstract — Sparse intersymbol-interference (ISI) channels are 
encountered in a variety of high-data-rate communication sys- 
tems. Such channels have a large channel memory length, but only 
a small number of significant channel coefficients. In this paper, 
trellis-based equalization of sparse ISI channels is revisited. Due 
to the large channel memory length, the complexity of maximum- 
likelihood detection, e.g., by means of the Viterbi algorithm (VA), 
is normally prohibitive. In the first part of the paper, a unified 
framework based on factor graphs is presented for complexity re- 
duction without loss of optimality. In this new context, two known 
reduced-complexity algorithms for sparse ISI channels are reca- 
pitulated: The multi-treUis VA (M-VA) and the parallel-trellis VA 
(P-VA). It is shown that the M-VA, although claimed, does not lead 
to a reduced computational complexity. The P-VA, on the other 
hand, leads to a significant complexity reduction, but can only be 
applied for a certain class of sparse channels. In the second part 
of the paper, a unified approach is investigated to tackle general 
sparse channels: It is shown that the use of a linear filter at the 
receiver renders the application of standard reduced-state trellis- 
based equaUzer algorithms feasible, without significant loss of op- 
timaUty. Numerical results verify the efficiency of the proposed 
receiver structure. 

Index Terms — Trellis-based equalization, sparse ISI channels, 
complexity reduction, prefiltering. 



1. Introduction 

SPARSE intersymbol-interference (ISI) channels are en- 
countered in a wide range of communication systems, such 
as high-data-rate mobile radio systems (especially in hilly ter- 
rain), wireline systems, or aeronautical/ satellite communica- 
tion systems. For mobile radio applications, fading channels 
are of particular interest. The equivalent discrete-time chan- 
nel impulse response (CIR) of a sparse ISI channel has a large 
channel memory length L, but only a small number of signifi- 
cant channel coefficients. 

Due to the large memory length, equaUzation of sparse ISI 
channels is a demanding task. The topics of linear equahzation 
and decision-feedback equalization for sparse ISI channels are, 
e.g., addressed in [1], where the sparse structure of the chan- 
nel is exphcitly utilized for the design of the corresponding 
finite-impulse-response (FIR) filter(s). TrelUs-based equaUza- 
tion for sparse channels is addressed in [2]-[4]. The complexity 
in terms of trellis states of an optimal trellis-based equalizer, 
based on the Viterbi algorithm (VA) [5] or the Bahl-Cocke- 
Jehnek-Raviv algorithm (BCJRA) [6], is normally prohibitive 
for sparse ISI channels, because it grows exponentially with 
the channel memory length . However, reduced-complexity 

^The VA is optimal in ttie sense of maximum-lilcelihood sequence estimation 
(MLSE) and ttie BCJRA in tlie sense of maximum a-posteriori (MAP) symbol- 
by-symbol estimation. Both algorithms operate on the same trellis diagram. All 
statements concerning complexity hold for both the VA and the BCJRA. 



algorithms can be derived by exploiting the sparseness of the 
channel. 

In [2], it is observed that given a sparse channel, there is 
only a comparably small number of possible branch metrics 
within each trellis segment. By avoiding to compute the same 
branch metric several times, the computational complexity is 
reduced significantly without loss of optimality. However, the 
complexity in terms of trelUs states remains the same and thus 
the storage expense. As an alternative, another equalizer con- 
cept coined multi-trellis Viterbi algorithm (M-VA) is proposed 
in [2] that is based on multiple parallel irregular trellises (i.e., 
time- variant trellises). The M-VA is claimed to have a signifi- 
cantly reduced computational complexity and storage expense 
without (much) loss of optimality. 

A particularly simple solution to reduce the complexity of 
the conventional VA without loss of optimality can be found 
in [3]: The parallel-trellis Viterbi algorithm (P-VA) is based on 
multiple parallel regular trellises. However, an application of 
the P-VA is only possible for a certain class of sparse channels 
having a so-called zero-pad structure. In order to tackle more 
general sparse channels with a CIR close to a zero-pad channel, 
it is proposed in [3] to exchange tentative decisions between 
the parallel trellises and thus cancel residual ISl. This modified 
version of the P-VA is, however, suboptimal and is denoted as 
sub -P- VA in the sequel. 

A generalization of the P-VA and the sub-P- VA can be found 
in [4], where corresponding algorithms based on the BCJRA 
are presented. These are in the sequel denoted as parallel-trellis 
BCJR algorithms (P-BCJRA and sub-P-BCJRA, respectively). 
Some interesting enhancements of the (sub-)P-BCJRA are also 
discussed in [4] . Specifically, it is shown that the performance 
of the sub-P-BCJRA can be improved by means of minimum- 
phase prefiltering [7], [8] at the receiver. A specific FIR ap- 
proximation of the infinite-length linear minimum-phase filter 
is used, which preserves the sparse structure of the channel. 
This guarantees that the sub-P-BCJRA can still be applied after 
the prefiltering. 

Alternatives to trellis-based equalization are the tree-based 
LISS algorithm [9] and the Joint Gaussian (JG) approach 
in [10]. In this paper, trellis-based equalization for sparse ISI 
channels is revisited. A unified framework based on factor 
graphs [11] is presented in Section II for complexity reduction 
without loss of optimality, and the M-VA as well as the P-VA 
are recapitulated in this new context. Specifically, it is shown 
that the M-VA does, in fact, not lead to a reduction of compu- 
tational complexity, compared to the conventional VA. More- 
over it is illustrated, why the optimal P-VA can only be appUed 
for zero-pad channels. In order to equalize general sparse ISI 
channels, a simple alternative to the sub-P- VA/ sub-P-BCJRA 



is investigated in Section III. For this purpose, the idea in [4] 
to employ prefiltering at the receiver is picked up. It is demon- 
strated that the use of a linear minimum-phase filter renders the 
apphcation of reduced-state equalizers such as [12], [13] fea- 
sible, without significant loss of optimaUty. The proposed re- 
ceiver structure is notably simple: The employed equalizer al- 
gorithms are standard, i.e., not specifically designed for sparse 
channel, because the sparse channel structure is normally lost 
after prefiltering. Solely the linear filter is adjusted to the cur- 
rent CIR (which is particularly favorable with regard to fading 
channels), where the filter coefficients can be computed accord- 
ing to standard techniques available in the hterature [7], [8]. In 
order to illustrate the efficiency of the proposed receiver struc- 
ture, numerical results are presented for various types of sparse 
ISI channels. Bit error rates (BERs) are achieved that deviate 
only 1-2 dB from the matched filter bound (at a BER of 10^'^). 
To the authors' best knowledge, similar performance studies for 
prefiltering in the case of sparse ISI channels have not yet been 
presented in the literature. 

II. Complexity Reduction without Loss of 
Optimality 

A general sparse ISI channel has a comparably large chan- 
nel memory length L, but only a small number of significant 
channel coefficients hg, g = 0, G <^ L, according to 

h := [ho 0^ hi 0^ ... 0^ hcf , (1) 

/o zeros /i zeros fc-i zeros 

with /i > integer for all i and L = X^j_o (/i -|- 1). A sparse 
ISI channel, for which /q = /i = ... = fc-i =■ f holds, is 
referred to as zero-pad channel [3]. 

Throughout this paper, complex baseband notation is used. 
The fc-th transmitted M-ary data symbol is denoted as x[k], 
where k is the time index. A hypothesis for x[k] is denoted 
by x[k] and a hard decision by x[k]. For simplicity, the channel 
coefficients are assumed to be constant over an entire block of 
data symbols (block length > L). The equivalent discrete- 
time channel model is given by 

G 

y[k] = hox[k] + ^hgx[k-dg] + n[k], (2) 

where y[k] denotes the fc-th received sample and n[k] the fc-th 
sample of a complex additive white Gaussian noise (AWGN) 
process with zero mean and variance a^. Moreover, 

a 

dg := J^^i-^ + 1)' 1 < 5 < G, (3) 

denotes the position of hg in h. 

In the sequel, the channel vector h is assumed to be known 
at the receiver. Moreover, an M-ary alphabet for the data sym- 
bols is assumed. The complexity in terms of trelUs states of 
the conventional Viterbi/ BCJR algorithm is given by 0{M^} 
and is thus normally prohibitive. Given a zero-pad channel, it is 
proposed in [3] to decompose the conventional trellis diagram 
with = M(-''+i)*^ states into (/ + 1) parallel regular trel- 
lises, each having only M*^ states. As will be shown in the 



sequel, a decomposition into multiple parallel regular trellises 
is not possible in the case of a more general sparse ISI channel. 

In order to decompose a given trellis diagram into multi- 
ple parallel trelhses, the following question is of central inter- 
est: Which symbol decisions x[k] are influenced by a certain 
symbol hypothesis .T[fco]? Suppose, a certain decision x[ki] is 
not influenced by the hypothesis i[fco], where fco and fci are 
two arbitrary (but flxed) time indices. Furthermore, let the set 
Xko {x[k] I x[k] depends on .T[fco]} contain all decisions in- 
fluenced by x[ko] and the set X^^ all decisions influenced by 
x[ki]. If these two sets are disjoint, i.e. fl X/.^ =0, the hy- 
potheses .T[fco] and x[ki] can be accommodated in separate trel- 
lis diagrams without loss of optimality. In this case, a decom- 
position of the overall trellis diagram into (at least two) paraUel 
regular treUises is possible. 

This fact is illustrated in the sequel for two example channels 
(L = 8 and G = 2 in both cases). The first channel is charac- 
terized by a CIR h(i) := [/lo /ii /i2]'^and 
flie second by a CIR h^^) := [/iq /ii /i2]'^. 
For the time being, a particular symbol hypothesis .x[fco] is con- 
sidered, and it is assumed that hard decisions x[k] are already 
available for all time indices k < ko. Moreover, it is assiuned 
that the hypothesis x[ko] does not have any impact on decisions 
x[k] with k > ko + DL {D > integer; in the example, we 
consider the case D = 2)^. Fig. 1 shows a diagram for the 
first example CIR. The diagram may be interpreted as a factor 
graph [11] and illustrates the dependencies between hypothe- 
sis x[ko] and aU decisions x[k], ko < k < ko + DL. It can be 
seen from (2) that only the received samples y[fco], j/[fco + 6], 
and y[fco + 8] are directly influenced by the data symbol x[ko]- 
Therefore, there is a direct dependency between the hypothesis 
x[ko] and the decisions x[ko], x[ko + 6], and i[fco + 8]. The re- 
ceived sample j/[fco+8], for example, is also directly influenced 
by the data symbol a;[fco + 2]. Correspondingly, there is as well 
a dependency between x[ko] and the decision a; [fco + 2]. The 
data symbols a;[fco+6] and a;[fco-|-8] again directly influence the 
received samples y[ko+l2], y[fco+14], and y[fco+16], and so on. 

Including all dependencies, one obtains the second graph of 
Fig. 1. As can be seen, there is a dependency between i;[fco] 
and all decisions x[ko + 2i'], where j/ = 0,l,..., [Z)L/2J. Con- 
sequently, in this example it is possible to decompose the con- 
ventional trellis diagram into two parallel regular trellises, one 
comprising only the time indices ko+2i' and the other one com- 
prising the time indices fco + 1 + 21^. While the conventional 
trellis diagram has trellis states, there are only M'' states 
in each of the two parallel treUises. (Moreover, a single trelUs 
segment in the parallel trellises spans two consecutive time in- 
dices.) This result is in accordance with [3], since the CIR h^^^ 
constitutes a zero-pad channel [/iq h[ /ij /13 h'^]^, 
where G = 4, / = 1, and h[ = /i2 = 0. Generally spoken, a 
decomposition of a conventional trellis diagram into multiple 
parallel regular trellises is possible, if all non-zero channel co- 
efficients of the sparse ISI channel are on a zero-pad grid with 
/ > 1. In this case, the optimal P-VA can be applied^; otherwise 
one has to resort to the sub-P-VA or to alternative solutions such 
as die M-VA. 

"^This corresponds to the assumption that a VA with a decision delay of DL 
symbol durations is optimal in the sense of MLSE. 

■^The P-VA is still optimal in the sense of MLSE. The finite decision delay 
DL is not required and has only been introduced for illustrative purposes. 
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Fig. 1. Dependencies between symbol hypothesis x[feo] and subsequent decisions £[fc]: Example CIR h(^' := [ho /ii h2]^ (L = 8, G = 2}. 



The second CIR h^^^ constitutes a counter example. Here, 
the symbol hypothesis x[ko] influences all decisions x[k'], 
ko<k' < ko+DL (not depicted due to space limitation). Con- 
sequently, a decomposition of the conventional trellis diagram 
into multiple parallel regular trellises is not possible here. Still, 
a decomposition into multiple parallel irregular trellises is pos- 
sible, as proposed in [2] in the context of the M-VA. However, 
a significant subset of the dependencies resulting from the cor- 
responding factor graph is neglected in [2]. If all dependencies 
are taken into account, the M-VA does not yield any complexity 
advantage over the conventional VA. The computational com- 
plexities in terms of the overall number of branch metrics that 
have to be computed for a single decision iffco] stated in Ta- 
ble 1, for the conventional VA, the P-VA (example CIR h^^'), 
and the M-VA (example CIR h^^)). 

III. Prefiltering for Sparse Channels 

The preceding section has shown that trellis-based equaliza- 
tion of a general sparse ISI channel is quite a demanding task: 
An application of the optimal P-VA (or the P-BCJRA) is only 
possible for zero-pad channels. In the case of a more general 
sparse channel, the suboptimal sub-P-VA with residual ISI can- 
cellation can be used. However, for a good performance the 
CIR should at least be close to a zero-pad structure [3]. The 
M-VA, on the other hand, was designed for sparse channels 
with a general structure, but does not offer any complexity ad- 
vantage over the conventional VA if all dependencies between 
the individual symbol hypotheses are taken into account. 

In order to tackle general sparse ISI channels, a simple alter- 
native to the sub-P-VA/ sub-P-BCJRA is proposed in the sequel: 
It is demonstrated that the use of a linear minimum-phase filter 
at the receiver renders the application of standard reduced-state 
equalizer algorithms feasible, without significant loss of opti- 
mality. The receiver structure under consideration is illustrated 
in Fig. 2, where z[k] denotes the fc-th received sample after pre- 
filtering and hmin the filtered CIR. Within the scope of this pa- 
per, the ideal linear minimum-phase filter is approximated by 



TABLE I 

Computational complexities: Viterbi algorithm (VA), 

PARALLEL-TRELLIS VA (P-VA), AND MULTI-TRELLIS VA (M-VA). 



Conventional VA, 

any CIR withL = 8 


P-VA, 

example CIR h*^' 


M-VA, 

example CIR h*^' 




0{2M^} 


0{3 M" + 2 M** 
+ 2M'^ + 2M'' 
+ 2 M'^ + 2 
+ 2M^ + 71/2} 



an FIR filter of length Lp, where the approach in [8] is used 
to calculate the filter coefficients. The resulting FIR filter ap- 
proximates a discrete-time whitened matched filter (WMF), i.e., 
the effect of noise coloring is negligible. The computational 
complexity of calculating the filter coefficients is 0{Lp L?), 
i.e. only linear with respect to the filter length. Therefore, 
comparably large filter lengths are feasible. With regard to the 
trellis-based equalizer, we focus on delayed decision-feedback 
sequence estimation (DDFSE) [13] in the sequel. The number 
of trellis states in the DDFSE equalizer is M^^ , where X <C i is 
a design parameter In order to obtain a complexity that is sim- 
ilar to that of the (sub-)P-VA/ P-BCJRA equalizer, one should 
choose K such that'* 



K < logM(/ + l) 



G. 



(4) 



In the following section, it is shown that the sparse channel 
structure is normally lost after prefiltering. Afterwards, numer- 
ical results are presented for various examples to demonstrate 
the efficiency of the proposed receiver structure. 

A. Channel Structure After Prefiltering 

The sparse structure of a given CIR h is normally lost after 
minimum-phase prefiltering. For example, consider the CIR 

"^In order to find an appropriate value for K in the case of a general sparse 
ISI channel, one has to iind an underlying zero-pad channel with a structure as 
close as possible to the CIR under consideration. 
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Fig. 2. Receiver structure under consideration. 

h=[\/o:5 Vol VQA]^ (i.e. I |h||2 = 1), and let i7(z) 
denote the z-transform of h. The zeros of H{z) are given by 

. zo 1.2 = 0.69 ± j0.80 (|zo i,2| = 1.06) and 

. zo 3,4 = -0.69 ± jO.56 (|zo 3,4! = 0.89). 
Correspondingly, h is neither minimum-phase nor maximum- 
phase^. The z-transform -ffmin(z) of the filtered CIR h,„in is 
obtained as follows: Those zeros of H{z) that are inside or on 
the unit circle are retained for iJmin(z), whereas the outside 
zeros are reflected into the unit circle. The resulting minimum- 
phase CIR is h„,in = [0.79 0.12 - 0.02 0.20 0.56]T, which 
is not sparse anymore. 

As opposed to this, in the special case of a zero-pad channel 
the sparse structure is always preserved after minimum-phase 
prefiltering: Let h. :— [h^ hi ... /ig ]^ denote a (non-sparse) 
CIR with z-transform H{z), and let hzp denote the correspond- 
ing CIR with memory length (/+1)G and z-transform Hzp{z), 
which results from inserting / zeros in between the individual 
coefficients of h. Furthermore, let zq, 1, zo. g denote the 
zeros of H{z). An insertion of / zeros in the time domain cor- 
responds to a transform z 1-^ z^/'^^'^^'^ in the z-domain, i.e., 
Hzp{z)= H{zf+^). This means, the (/+1)G zeros ofiJzp(z) 
are given by the (/+!) complex roots of zq, 1, zq. g. respec- 
tively. Consider a certain zero zij g :— t'o, g exp(j i/Jo. g) of 
H{z) that is outside the unit circle (ro, g > 1). This zero will 
lead to (/+!) zeros 

4'1 --^fr^exp(j?!^^±|M) (5) 

of Hzp{z) (A = 0,...,/) that are located on a circle of ra- 
dius T'o^g'^^^^ ^ 1' ^1^° outside the unit circle. By means 
of minimum-phase prefiltering, these zeros are reflected into 
the unit circle, i.e., the corresponding zeros of -ffzp.minl^) are 
given by I/zq'^^*. Therefore, the sparse CIR structure is re- 
tained after minimum-phase prefiltering (with the same zero- 
pad grid), since the zeros of iJzp,min(2^) are the (/+!) roots 
of the zeros of H^^iniz). Specifically, the non-zero coefficients 
of hzp,min are given by the CIR hnii„. If the zeros of H{z) (or 
equivalently of i?zp,min(2)) are not too close to the unit circle, 
hmin is characterized by a significant energy concentration in 
the first channel coefficients. In this case the effective channel 
memory length of hzp is significantly reduced by minimum- 
phase prefiltering, namely by some multiples of (/+!), cf. (1). 

B. Numerical Results 

In the sequel, numerical results obtained by Monte-Carlo 
simulations are presented. To start with, a static sparse ISI 
channel is considered, and the BER performance of the pro- 
posed receiver structure is compared with the sub-P-BCJRA 
equalizer [4]. As an example, we consider here the CIR 

^Given a minimum-pliase (maximum-phase) CIR h, all zeros of H{z) are 
inside (outside) the unit circle, i.e., I^o.i I < 1 (|2o,i | > 1) for all i. 
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Fig. 3. BER performance of the proposed receiver structure in the case of a 
static sparse ISI channel. 

h = [ /iQ . . . /i4 . . . /i7 . . . /ii5 with ha = 0.87 and 
/i4 = hj ^ /ii5 = 0.29 from [4], which has a general sparse 
structure (i.e., no zero-pad structure). The BER performance 
(binary antipodal transmission, M = 2) of the sub-P-BCJRA 
equalizer and the DDFSE equalizer with WMF is displayed in 
Fig. 3, as a function of E'b/A^o in dB, where i?b denotes the 
average energy per bit and A^o the single-sided noise power 
density {E^/Nq := l/Cn). Due to the given channel memory 
length, the complexity of MLSE detection is prohibitive. As a 
reference curve, however, the matched filter bound (MFB) is in- 
cluded, which constitutes a lower bound on the BER of MLSE 
detection. The filter length of the WMF has been chosen as 
L-p = 40. Since the channel is static, the filter has to be com- 
puted only once. When the parameter K is chosen as if = 4, 
the overall receiver complexity is approximately the same as for 
the sub-P-BCJRA equalizer As can be seen in Fig. 3, the BER 
performance achieved by the proposed receiver structure is also 
comparable to the sub-P-BCJRA equaHzer. At a BER of 10^^, 
the loss with respect to the MFB is only about 1 dB. At the ex- 
pense of a small loss due to residual ISI (0.5 dB at the same 
BER), the complexity of the DDFSE equalizer can be further 
reduced to if = 3. 

Next, we consider the case of a sparse Rayleigh fading chan- 
nel model, i.e., the channel coefficients hg (g = 0, ...,G) in (1) 
are now zero-mean complex Gaussian random variables with 
variance E{|ft,gp} =: cr^ . It is assumed in the sequel that 
the individual channel coefficients are statistically independent. 
Moreover, block fading is considered for simplicity. As an ex- 
ample, we assume a CIR with G = 3 and a power profile 

P := [<o 0_^ <i crl^ al,]^. (6) 

/ zeros 

Note that this CIR again does not have a zero-pad structure. By 
choosing different values for the parameter /, different chan- 
nel memory lengths L = f + 6 can be studied. To start with, 
consider a power profile with equal variances cr^ q = ... = 
cr^ 3 = 0.25. Fig. 4 shows the BER performance of the pro- 
posed receiver structure for binary transmission and three dif- 
ferent channel memory lengths L (solid lines: L = 6, dashed 
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Fig. 4. BER performance of the proposed receiver structure in the case of a 
sparse Rayleigh fading channel. 

lines: L = 12, dotted lines: L — 20). The filter length has been 
chosen as Lp = 20 (i = 6), Lp = 36 (i = 12), and ip 60 
(L = 20). As reference curves, the BER for flat Rayleigh fading 
(L = 0) is included as well as the MFB [14, Ch. 14.5]. In the 
case L — Q, MLSE detection is still feasible. As can be seen in 
Fig. 4, its performance is very close to the MFB. The DDFSE 
equalizer with K = 5 in conjunction with the WMF achieves 
a BER performance very close to MLSE detection (the loss at 
a BER of 10^"^ is only about 0.6 dB). Even when the chan- 
nel memory length is increased to i = 20, the BER curve of the 
DDFSE equalizer with WMF deviates only 2 dB from the MFB. 
Specifically, a significant gain compared to flat Rayleigh fading 
is achieved, i.e., a good portion of the inherent diversity (due 
to the independently fading channel coefficients) is captured. 
When the DDFSE equalizer is used without WMF, a significant 
performance loss occurs already for L = 6. For the larger chan- 
nel memory lengths, the influence of the WMF makes a dra- 
matic difference: The BER increases by several decades when 
the WMF is not used. 

Similar performance results have also been obtained for un- 
equal variances af-^ g. When the power profile p of the original 
CIR already exhibits an energy concentration in the first channel 
coefficients, the benefit of the WMF is smaller, but still signif- 
icant. As an alternative to the proposed receiver structure, we 
have also studied the use of a linear channel shortening filter 
(CSF) [15], which transforms a given CIR with memory length 
L into a CIR with a reduced memory length Ls<L. When Ls 
is chosen small enough, the application of the conventional VA 
is again feasible, operating on the shortened CIR. However, the 
performance of this receiver structure turned out to be inferior 
to the DDFSE equaUzer with WMF [16]. 

The concept of minimum-phase prefiltering for sparse ISl 
channels is also beneficial when using a tree-based equalization 
algorithm, such as the LISS algorithm [9]. In order to obtain a 
small complexity, the metrics of two competing paths that devi- 
ate closely to the root of the tree should differ as much as possi- 
ble. This is achieved by means of minimum-phase prefiltering, 
due to the energy concentration in the first channel coefficients 
of the filtered CIR. 



IV. Conclusions 
In this paper, trellis-based equalization of sparse inter- 
symbol-interference channels has been revisited. Due to the 
large memory length of sparse channels, efficient equalization 
with an acceptable complexity-performance trade-off is a de- 
manding task. With regard to complexity reduction, it has been 
demonstrated in which cases a decomposition of the conven- 
tional trellis diagram into multiple parallel trellises is possible 
without loss of optimality. In order to tackle general sparse 
channels, a receiver structure with a linear filter and a reduced- 
complexity trelUs-based equalizer has been studied. The em- 
ployed equalizer algorithm is standard, i.e., not specifically de- 
signed for sparse channels, because the sparse channel structure 
is normally lost after prefiltering. The coefficients of the linear 
filter can be computed using standard techniques from the liter- 
ature. By means of numerical results, the efficiency of the pro- 
posed receiver structure has been demonstrated, both for static 
and fading channels. 
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